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^ ABSTRACT 

We combine variability information from the MAssive Compact Halo Objects (MA- 
CHO) survey of the Large Magellanic Cloud (LMC) with infrared photometry from the 
Spitzer Space Telescope Surveying the Agents of a Galaxy's Evolution (SAGE) survey 
to create a dataset of ^■SOOOO variable red sources. We photometrically classify these 
sources as being on the first ascent of the Red Giant Branch (RGB), or as being in 
one of three stages along the Asymptotic Giant Branch (AGB): oxygen- rich, carbon- 
rich, or highly reddened with indeterminate chemistry ("extreme" AGB candidates). 
We present linear period-luminosity relationships for these sources using 8 separate in- 
frared bands (J, H, Kg, 3.6, 4.5, 5.8, 8.0, and 24 nm) as proxies for the luminosity. We 
find that the wavelength dependence of the slope of the period-luminosity relationship is 



different for different photometrically determined classes of AGB stars. Stars photomet- 
rically classified as 0-rich show the least variation of slope with wavelength, while dust 
enshrouded extreme AGB stars show a pronounced trend toward steeper slopes with 
increasing wavelength. We find that 0-rich AGB stars pulsating in the fundamental 
mode obey a period-magnitude relation with a slope of —3.41 it 0.04 when magnitude 
is measured in the 3.6 fiui band, in contrast to C-rich AGB stars, which obey a relation 
of slope -3.77 ±0.05. 

Subject headings: (galaxies:) Magellanic Clouds, infrared: stars, stars: AGB, stars: 
carbon, stars: variables: general 



1. INTRODUCTION 

As intermediate mass stars (~l-8 M©) exhaust the helium in their cores during the late stages 
of stellar evolution, they begin their ascent of the Asymptotic Giant Branch (AGB). These stars 
consist of an inert C/0 core surrounded by concentric shells of helium and hydrogen. During 
the early AGB (E-AGB) phase, the evolution is driven by shell hydrogen burning. Stars are 
typically oxygen-rich during the E-AGB phase, with a photospheric C/0 ratio < 1. Eventually, the 
helium shell ignites, beginning the thermally-pulsing (TP-AGB) phase. Subsequent evolution on 
the AGB consists of relatively long periods of hydrogen burning, punctuated at regular intervals by 
brief helium shell flashes (thermal pulses) that rapidly convert the helium to carbon, dramatically 



increasing the star's brightness for a brief period (Vassiliadis & Wood 1993 Habing & Olofsson, 
editors, ||2003 >. These thermal pulses drive large-scale convective zones in the stellar interior which 
can "dredge-up" nuclear-processed material, most notably newly created carbon, to the stellar 
surface, leading to changes in the observable surface chemistry of the star. This is the third dredge- 
up for intermediate mass stars, and the result is that for stars with M < 4 Mq, the surface C/0 
ratio changes to become > 1 at which time the star is referred to as carbon-rich. For stars with 
M > 4 M0, the temperature at the base of the convective shell is sufficiently high (~ 10^ K) to 
burn the C into N, a process known as "hot-bottom burning (HBB)," and the star remains 0-rich. 

In addition, hydrodynamic instabilities drive shock waves through the stellar interior, levi- 
tating the outer layers of the star into cooler regions where dust grains can condense. Radiation 
pressure then drives the dust (and the gas to which it is collisionally coupled) into the interstellar 
medium (ISM). These pulsations cause AGB stas to exhibit complicated, multi-periodic variations 



in brightness (Eraser et al. 2005). In addition, this process can cause an AGB star to lose mass at 



rates up to ~ 10 ^ Mq yr ^ (van Loon et al. 1999). This makes AGB stars one of the dominant 



sources of dust in the universe, and a key component in the chemical evolution of galaxies. Toward 
the end of its AGB lifetime, this mass loss can completely enshroud a star with a thick circumstellar 
dust shell. These "extreme" AGB stars are nearly undetectable in the optical, and because the light 
from the star is dominated by thermal emission from the circumstellar shell, cannot be classified 



as either 0-rich or C-rich based on their near-IR photometry (Cioni et al. 2006). Most of them 



are beheved to be C-rich, but the brightest should be so massive that HBB has left them 0-rich 



(Matsuura et al. 2009). Without spectroscopic confirmation, we do not propose classifications for 



specific sources of this class in this paper. 



The Surveying the Agents of a Galaxy's Evolution (SAGE) survey (Meixner et al. 2006) is an 
unbiased 7° x 7° survey of the Large Magellanic Cloud (LMC) using the IRAC (3.6, 4.5, 5.8, and 



8.0 //m) and MIPS (24, 70, and 120 //m), instruments aboard the Spitzer Space Telescope (Werner 
et al. [2004 ), intended to trace the life-cycle of the baryonic matter in the LMC. This wavelength 
regime makes SAGE extremely well suited to studies of evolved stars. The LMC makes an ideal 
target for studies of stellar populations, as its high galactic latitude minimizes the foreground 
contamination, and its distance of ~50 kpc (van Leeuwen et al. 2007) makes it possible to resolve 



individual stars while simultaneously being able to neglect its 3-dimensional structure. 

Mira variables (a subset of AGB stars) have been known to follow a linear relationship in 
magnitude, log Period space for nearly a century ( Gerasimovic|1928 Feast et al. 1989, and references 



therein). Wood & Sebo (1996) extended the known dataset of AGB variables to shorter periods and 
established that Mira-type variables pulsate in the fundamental mode, and that other, higher order 
pulsational modes exist as well. Expanding the number of data points to ~1000, Wood et al. ( 1999 ) 



identified 4 parallel linear sequences in period-luminosity (P-L) space, and proposed mechanisms 
for each. Ita et al. (2004), using a sample of ~30 000 variable evolved stars with period information 



derived from the Optical Gravitational Lensing Survey (OGLE;|Udalski et al.|[l997), further refined 



these sequences (see section 2.1). Also using OGLE data, Soszynski et al. (2004) identified a new 



P-L sequence in the secondary, non-dominant variation in the light curves of evolved stars, and 
proposed this as a means of separating stars dimmer than the tip of the red giant branch (TRGB) 
into stars on the first ascent of the RGB and those evolving along the AGB. With a dataset of 



comparable size derived from the MACHO survey (Alcock et al. 1997), Fraser et al. (2005, 2008) 



further investigated the multi-periodic nature of stars on the sequences identified by Wood et al. 



(1999). Using early results from the SAGE survey. Glass et al. (2009) extended work on the P-L 



relation in AGB stars by investigating the dependence of the calculated slope of the relation on 
wavelength. 

Our survey combines the SAGE project's photometry archive and the MACHO variability 
dataset to produce the largest sample of variable evolved stars, with photometry measured in 
multiple near-IR wavelengths, to date. The unprecedented size of our dataset, combined with our 
extensive spectral coverage, reveals new patterns and details in the relationship between luminosity, 
variability, and spectral energy distribution amongst evolved stars. 

This paper is organized as follows: § [2] details our source selection process, the identification 
of the multiple sequences observed in period-luminosity space, and the division of our sample 
into 6 photometrically determined categories. §[3] presents the quantitative and qualitative effects 
observed within our sample. § |4] compares our current study to previous work in the field, and § [5] 



summarizes the conclusions of the current work. 



2. DATA 



The SAGE survey observed the LMC in two epochs, three months apart. From each epoch, 
two source hsts, an archive and a catalog, were generated. The catalog stresses reliability over com- 
pleteness and requires more stringent data quality standards (higher signal to noise, less crowding 
in the field) for inclusion of a source than the archive. The standards for source inclusion in both 



archive and catalog are detailed in Meixner et al. (2006). 



In 2009, the SAGE team released the SAGE mosaic photometry catalog and archive to the 
Spitzer Science Center. By combining data from all epochs of the SAGE survey, the mosaic pho- 
tometry catalog and archive are more complete than the individual epochs, and photometric errors 
are reduced. The mosaic photometry data are well documented in the SAGE Data Products De- 
scription documenlj^ The SAGE Mosaic Photometry Catalog contains ~6.5 million sources, and 
the archive contains ~7 million. 



We extract a list of 32 744 sources from Eraser et al. (2008)'s catalog of evolved Long Period 



Variables (LPVs) in the LMC with a counterpart in the SAGE Mosaic Photometry Archive. We 
utilize the archive in order to maximize sample size. We use a matching radius of 2", and keep only 
the closest SAGE source to a given MACHO source. We do not include the ~11 000 sources [Eraser 



et al. (2008) identify as the "one- year artifact," which have an artificial period of 365 days caused by 



MACHO 's annual observing schedule. The database merging process can create duplicate entries 
for some sources, due to slight (< 0.1 ") shifts in position between SAGE Epoch 1, Epoch 2 and the 
single frame mosaic photometry. By requiring that every source in our list have a unique MACHO 
field. tile. sequence identifier, we find ~500 of these duplicate matches in our dataset and cull them 



from our final list. Based on the definitions used by Eraser et al. (2008), we assign stars to one of 6 



roughly parallel sequences in period- i^Tg space (Eig.[T]). Our definitions are detailed in Appendix [A} 

We are left with a final dataset of 30 747 evolved stars with well-defined MACHO periods 
and high-quality near IR photometry from 2MASS and SAGE. Table [l] compares the number of 
sources in the present study to those in Eraser et al.| ( |2008[ ), Srinivasan et al. (2009), and Vijh et 
al. (2009). Our entire dataset is available online. Table ^ presents the photometric information for 
a few sources from our sample as a guide to the format of the online table. 



L 



http : //data, spitzer . caltech . edu/populcir/sage/20090922_enhanced/documents/ 
SAGEDataProductsDescription_Sep09 .pdf 



2.1. Period-Luminosity Sequences 



Wood et al. ( 1999 ) identified 5 parallel period-luminosity sequences in the MACHO dataset and 



proposed underlying physical mechanisms for them. In order of increasing period, the sequences of 
pulsating stars were named A, B, and C, with sequences D and E exhibiting variation due to an 
unknown mechanism. Sequence C was identified as pulsation in the fundamental mode, sequence 



B as the first and second overtones, and sequence A as the third overtone. Kiss Sz Bedding (2003) 



and Ita et al. (2004) split sequence B into two sequences (C and B) representing pulsation in 



the first and second overtone respectively. Eraser et al. ( 2005 ) retained the names of sequences D 



and E from Wood et al. (1999), but renamed the pulsation sequences 1-4 in order of decreasing 



period, in order to easily accommodate shorter period sequences, and align the naming convention 



with increasing pulsation overtone. We follow the naming convention of Eraser et al. (2005, 2008) 



because of this parallel with theory. Numbering the sequences in order of decreasing period means 
that higher sequence numbers generally correspond to higher order pulsational modes, aligning the 



empirical naming convention and a theoretical explanation more gracefully. Nicholls et al. (2010) 



has provided convincing evidence that sequence E consists of ellipsoidal binaries, systems in which 



the red giant member has filled its Roche lobe. Recent work (e.g Nie et al. 2010, Nicholls et al. 



2009) has demonstrated that the variability seen on sequence D cannot be due to binarity, but 



its actual cause remains unknown. We keep the alphabetic names for these sequences to separate 
them from the stellar pulsation sequences. Table [3] summarizes these conventions, and Eigure [l] 
illustrates these sequences. The precise definitions of these sequence are discussed in Appendix \K\ 



2.2. Identification and Classification of AGE Candidates 



Before being photometrically classified, all sources were de-reddened in the J, H and Kg bands 
following Glass (1999). Specifically, we use corrections of Aj = 0.112, Ah = 0.065, and Ak = 



0.037 mag. We define AGB candidates using 2MASS and SAGE photometry (Cioni et al. 2006 



Blum et al. 2006: Srinivasan et al. 2009). We identify a source as an AGB candidate by requiring 



that 

Ks > -13.333 X (J - Ks) + 24.666 and Ks < 12.05 (1) 

We classify sources as low- or moderately-obscured oxygen-rich and carbon-rich AGB candidates 
based on their near-IR photometry. Specifically, a source is classified as an 0-rich candidate if it 
lies leftward (blueward) of the line 



Ks 



-13.333 X {J -Kg) + 28.4: 



(2) 



in the Kg vs. J — Kg color-magnitude diagram (CMD), and a C-rich candidate if it lies rightward 
(redward) of this line ( Cioni et al.|2006 fig. 1). Extreme AGB stars are defined using their position 
in the [3.6] vs. J- [3.6] CMD 



J- [3.6] > 3.1 and [3.6] < 10.5 



(3) 



or the IRAC [8.0] vs. [8.0] — [3.6] CMD if a source lacks a valid J band magnitude. These sources 
are thought to be highly evolved AGB stars, enshrouded by thick circumstellar dust shells with 
high rates of mass loss. 



The TRGB is located at K^ = 12.0 (e.g. Srinivasan et al. 2009) or / = 14.54 (Cionietal 



|2006D, and is typically identified by a clear fall-off in population density (e.g. Nikolaev &: Weinberg 
|2000p . However, while the TRGB is justified on physical grounds as a maximum luminosity for 
RGB stars, it does not represent a physically motivated minimum luminosity for stars on the AGB. 



Examination of a plot of Kg magnitude vs. variability period ( Fraser et al.[2008 fig. 2), shows that 
the sudden decrease of population density brighter than K^ = 12.0 indicative of the TRGB only 
appears on Period-Luminosity sequences 3 and 4. We thus classify a star as an RGB candidate 
if it is fainter than the K^ or / band TRGB cutoff, and it lies on either sequence 3 or 4. Stars 



dimmer than the TRGB on sequences 1 and 2 are classified as 0-rich AGB candidates. Soszynski 



et al. (2004) proposed a means of separating AGB stars dimmer than the TRGB on sequence 3 



or 4 from the dominant RGB population on these sequences. Their method relied on analysis 
of non-dominant periods and a distinct P-L sequence with shorter periods than sequence 4 (this 
sequence would be called sequence 5 in our nomenclature). However, we do not detect this sequence 
in either the dominant or secondary pulsation periods of our sample, and thus we do not employ 
this method. The use of periods in addition to the dominant two periods available to us grants 



Soszynski et al. ( 2004 ) finer resolution in separating RGB stars from AGB candidates dimmer than 



the TRGB, but the consistency of our results with those of other authors gives us confidence in 
the broad correctness of our classification. Figure [2] reproduces Figure [T| but this time color coded 
according to the photometrically determined chemical classification we assign to the stars in our 
sample. Table [4] details the number of each class of star on each P-L sequence. 

Figure [SJ a histogram of [3.6] magnitude color-coded by P-L sequence, provides further justifi- 
cation for our classification. Sequences 3 and 4 are dominated by a large population dimmer than 
[3.6] w 12. Since there is no theoretical expectation for AGB stars to pile up at the TRGB, we 
interpret this as evidence for a large RGB population on Sequences 3 and 4. Sequences 1 and 2 
do show evidence of a distinct population of stars at this brightness (more pronounced in sequence 
2), but neither shows the dramatic RGB population found in sequences 3 and 4. We interpret this 
as a small amount of RGB contamination in sequences 1 & 2. All four sequences show a stellar 
population with a peak at [3.6] ~ 11.2, coincident with our photometrically determined 0-rich 
AGB candidate population. In addition, the stars photometrically classified as C-rich stars are 
visible in sequences 1 and 2 as a peak in the magnitude distribution brighter than [3.6] ~ 10. This 
population is also visible in sequence 3 but absent from sequence 4. These trends in the luminosity 
functions of the P-L sequences mirror and lend support to the photometrically determined chemical 
classifications detailed in Table HI 



Highly evolved AGB stars and young stellar objects (YSOs) have similar mid-IR colors. Whit- 

it the region of the [8.0] vs. [8.0] -[24] CMD defin 

[24] > 2.2 and [8.0] > 11 - 1.33 x ([8.0] - [24]) 



ney et al.| f2008| ) suggest that the region of the [8.0] vs. [8.0] -[24] CMD defined by 

[8.0] 



constitutes a region in color-magnitude space dominated by YSOs. As all of our sources have well 
determined MACHO periods, none of our sample are YSOs. Of our ~30 000 stars, we find that 
only 150 sources (109 RGB, 3 extreme AGB and ~40 0-rich and C-rich AGE stars) fall in this 
region. We note that the optical wavelengths used by the MACHO survey biases our dataset to 
under-sample this region of CMD space, but the extremely small fraction of our sample that falls 



into this region lends support to the criteria proposed by Whitney et al. (2008) 



Our sample consists of 17 059 AGB candidate stars: 12172 sources photometrically classified 
as 0-rich, 4455 as C-rich, and 432 as extreme AGB candidates. Based on near-IR photometry and 
variability period, 13 688 sources are classified as stars on the RGB. 



The Ks magnitude is often used as a proxy for luminosity for AGB stars (e.g. Wood et al.|1999 



Cioni et al. 2006). The inclusion of Spitzer IRAC and MIPS data at wavelengths out to 24 //m 



reveals new features in the period magnitude diagram, especially among the reddest sources, the 
extreme AGBs, whose SEDs peak at 3 jjlvo. or redward ( Vijh et al.|2009 ). FigH^compares the 2MASS 
Ks band, the IRAC 4.5, 5.8, and 8.0 //m bands, along with the MIPS 24 /xm band as luminosity 
proxies, all plotted against period as in Figs. [l]&:[2] 

The Ks band, in the upper left, is a reasonable approximation for the brightness of the low- 
obscuration RGB, 0-rich and C-rich candidates, but comparison with the other panels in Fig. |4j 
shows that this approximation fails for the extreme AGB candidates. The RGB, 0-rich and C-rich 
candidates appear in roughly the same positions, relative to to one another, in the 2MASS and 
IRAC bands. The extreme sources, however, clearly have their luminosity underestimated by the 
Ks band. These sources, the most luminous in our sample ( Srinivasan et ali]|2009 ) , appear no more 
luminous than the C-rich candidates at this wavelength. The IRAC bands, on the other hand, 
correctly place the extreme candidates at the brightest end of sequence 1, the fundamental mode 
pulsators. 

The 24 //m band is sensitive to emission from cool dusty envelopes. The observed flux from 
the sources surrounded by the most optically thick dust shells is dominated by the emission from 
their outermost, coolest layers. For this reason, they are able to produce significant 24 ^m fiuxes 
(see Fig. 10, Srinivasan et al. 2009). The RGB/AGB stars with optically thin dust shells are 
either faint or undetected at 24 //m (for example, < 3% of the RGB stars in our sample have valid 
measured 24 fiva fiuxes) because their SEDs peak at much shorter wavelengths. The 24 ^m fiux 
is therefore not representative of the intrinsic brightness of these sources. This poor estimate for 
intrinsic brightness causes the sequence structure to break down entirely. 



3. RESULTS 
3.1. Period Magnitude Relations 

Using the method of weighted least-squares, linear models of the form 

m = AxlogP + B (4) 

were fit to each of the sequences as a whole, and to each subgroup (0-rich, etc.) within each 
sequence. Periods and amplitudes are from the dominant mode of variability derived by [Eraser 



et al. (2005, 2008) from MACHO blue-band lightcurves. Periods were considered free of errors, 
and because the SAGE observations are not phase corrected, magnitudes were weighted by the 
inverse of the quadrature sum of the 1-cr photometric error and one-half the MACHO peak-to-peak 
variability amplitude. Table [5] lists the linear fit parameters thus determined for each of the 6 
sequences, using the IRAC [3.6] magnitude as the luminosity proxy. 

Previous work has indicated that Mira variables exhibit much lower amplitude variation in 



the IR than in the optical (Lattanzio & Wood (2003), in: Habing & Olofsson, editors, 2003), and 
thus using one-half of the MACHO amplitude in the error term might overestimate the true scatter 
in the data due to intrinsic variation. To investigate this, we extracted a sample of ~100 sources 



from the catalog of Vijh et al. (2009) with well-determined MACHO periods between 150 and 200 
days. At these periods, the three month cadence of the SAGE observations could sample these stars 
anywhere between maximum and minimum brightness. Thus the range of variation in brightness 
between Epoch 1 and 2 observations could represent the full mid-IR peak to peak amplitude of 
variation. The MACHO amplitudes for the 0-rich and C-rich AGB candidates are statistically very 



similar to the variation seen in this smaller sample. The small sample from Vijh et al. (2009) has 
a median mid-IR variation of ~0.2 mag between SAGE Epoch 1 and Epoch 2. The median optical 
amplitudes for all 0-rich and C-rich sources together are 0.2 and 0.3 magnitudes respectively. In 
addition, we compared the distribution of the MACHO variability amplitudes to the distribution 
of the residual scatter (Table ^ about the P-L fits. The residual scatter of ~0.5 magnitude about 
the best fit line (see Table ^ is due to the intrinsic variability of the sources. The SAGE mosaic 
photometry is constructed from the two epochs of the SAGE survey, but lacks the time-stamp 
information necessary for the phase-correction necessary to compute a true mean magnitude for 
these variable sources. Some of our sources have been observed near maximum brightness, while 
other were near their minimum brightness. The residual scatter about the P-L relations gives us 
an indication of the scale of the IR variability of our sample. Both quantities followed similar 
distributions, and were of the same scale. This justifies our use of the MACHO optical variability 
amplitude as a measure of the IR variability of our sources for the estimation of the uncertainty in 
the IR mean magnitudes in the P-L relations. 

Figureplshows the linear fits to the stars in sequence 1 in the [3.6] vs. period diagram. Sequence 
1 is consistent with stellar pulsation in the fundamental mode (Eraser et al.||2005), and was chosen 



because it is the only sequence which contains a significant number of detected extreme AGB 



candidates. The 3.6 //in IRAC band was chosen because it is the most sensitive of the SAGE bands 
( Meixner et al.|2006 ) . The different trends followed by the different types of evolved star candidates 



is apparent. The plotted linear fits are: 

0-rich [3.6] = (-3.41 ± 0.04) x log P + (18.88 ± 0.09) 

C-rich [3.6] = (-3.76 ± 0.05) x log P + (19.35 ± 0.12) 

Extreme AGB [3.6] = (-4.27 ± 0.19) x log P + (20.4 ± 0.49) 

The increasing steepness of the linear relations followed by the different categories of AGB star 
would be consistent with a non- linear P-L relation if the evolved variable stars were not divided into 
sub-categories. We experimented with fitting quadratic functions to all the stars in each sequence, 
regardless of chemical classification. While we found such fits to be consistent with our data, we 
did not find them to be statistically superior fits to the linear models presented here. 

Table [5] lists the parameters of the linear fits to each of the 6 sequences in the 3.6 fim band. 
Identical tables for all of the 2MASS, IRAC and MIPS bands used in this work are included as 
Appendix [B| in the online edition. 

Figure |6] plots the slope of the linear fit to the sources of each type of evolved star, in sequences 
1 and 2, in all the IR bands. Because we only use the period measured from the light curve in the 
MACHO blue band, this is a fixed value for all of our sources. The slope of the log P-magnitude 
relation is therefore a measure of how different the brightness in a particular band is between 
the longest period stars of a particular type and the shorter period stars of the same type on 
the same sequence. Figure [6] shows that 0-rich AGB candidates follow roughly the same linear 
relationship in all IR bands, showing that there is little difference between the SEDs of the 0-rich 
AGB candidates with the longest periods, and those with the shortest periods. The extreme AGB 
candidates, on the other hand, show a very strong trend toward steeper slopes as one moves towards 
longer wavelengths. The slope of the P-L relation is a measure of the dependence of brightness in a 
certain band to period. Smaller values for slope indicate smaller differences in brightness between 
short and long period stars within one class on one sequence (e.g. RGB stars on sequence 3). 
Steep slopes indicate large differences between short and long period stars within one class on one 
sequence (e.g. extreme AGB stars on sequence 1). All extreme AGB candidates become brighter 
as one looks at redder bands because they possess dusty circumstellar envelopes which reprocess 
the near-IR stellar photospheric emission to longer IR radiation. However, this mid-IR brightening 
is more dramatic for the longest period sources on a given P-L sequence, which causes the slope to 
become steeper and indicates that the longest period sources have the thickest circumstellar dust 
shells. Thus, the steep slopes of the P-L relations followed by extreme AGB candidates at longer 
wavelengths indicate not just the presence of circumstellar matter, but show that period is more 
closely tied to the presence of circumstellar matter in extreme AGB stars than in, e.g. 0-rich AGB 
stars. 

The slope of the log P-magnitude relation for C-rich AGB stars shows a more complicated 
variation with wavelength. In the IRAC [5.8] band (and to a lesser extent in the 4.5 /im band). 
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the slope becomes noticeably less steep, indicating that there is less brightness contrast between 
the C-rich stars with longest periods and those with shorter periods on the same sequence in this 



band. The spectra of C-rich stars have a strong CO absorption feature at ~4-5 //m (e.g. Aringer et 



al. 2009 fig. 4). The observed decrease in the slope of the log P-magnitude relation in the 5.8 //m 
band compared to neighboring bands indicates a possible relationship between the strength of this 
feature and the pulsation period of carbon-rich stars. Specifically, that longer-period C-rich AGB 
stars may have stronger CO absorption features than those stars with shorter periods. While many 
of them are thought to be carbon-rich, extreme AGB candidates do not show this effect because 
their IR emission is dominated by their circumstellar dust shells. The absence of this effect in O- 
rich AGB candidates, which also feature CO in their photospheres, argues against our hypothesis. 
Further IR spectroscopic follow-up, such as the SAGE-SPEC survey ( Kemper et al.|20"l0 ) will shed 
further light on this. 

Figure |6] also has implications for future observational strategies. AGB stars are a dominant 
source of IR light in galaxies with intermediate- and old-age stellar populations. With major IR 
observatories such as the James Webb Space Telescope coming online in the next few years, the AGB 
period- luminosity relationship constitutes a possible distance indicator. Figure [6] shows that the 
observed scatter in the P-L relation could be reduced, while maximizing sample size, by observing 
in the near-IR, specifically the Kg band (~2.2 iim). In this band, all AGB candidates brighter than 
the TRGB obey very nearly the same P-L relation. Furthermore, the slope of the relation measured 
for AGB candidates in the Kg band will be the slope of the relation followed by O-rich AGB stars 
(which can be photometrically classified) at all IR wavelengths out to 8 /im (Also see Table p| . 
Figure [8] plots the median SED for each classification of star used in this paper. Error bars represent 
the spread of each flux within the population, not measurement uncertainty. The Kg band is the 



location of the peak of the SED for the majority of C-rich AGB stars (also see Srinivasan et al 



2009), and is very near the peak for both the O-rich AGB stars and those on the first ascent of 
the RGB. Furthermore, 97% of our sample (including 82% of the extreme AGB stars) have a valid 
flux measured in the Kg band, implying that this band will provide an adequate sample size. The 
intrinsic brightness of the reddest, most evolved AGB stars will be severely underestimated (see 



2.2), but these stars can be identified by near-IR color and excluded from the sample to reduce 
scatter in the derived relationship. Table |6] presents all the derived Kg vs. logP relations for our 
sample. 



3.2. Long Secondary Periods 

Approximately 30% of variable AGB stars exhibit long secondary periods (LSPs), variations 
which occur on time scales nearly an order of magnitude longer than pulsation in the fundamental 
mode (Sequence 1). The LSP phenomenon is represented in our data set by sequence D, and the 



mechanism behind it is still unknown (Nicholls et al. 2009 Nie et al. 2010). Eraser et al. (2008) 



noted that the variability properties of sequence D are consistent with a population drawn from all 
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of the other P-L sequences. The LSP is known to be related to mass- loss ( Wood &: Nicholls||2009 ) , 
but in an as yet unknown manner. 

The lower left panel of figure [7] shows the variation of the LSP with wavelength. The slopes of 
the P-L relations followed by the AGB candidate stars on sequence D show the same dependence 
on wavelength as the stars on the other sequences. 0-rich AGB stars show little variation in P-L 
relation slope with wavelength, and C-rich AGB stars show identical behavior in Sequences 1, 2, 3 
and D. We do not include the Extreme AGB candidates from sequence D in Figure [7] because the 
fits, based on only 18 stars, are not statistically significant (Table rol). Whatever mechanism lies 
behind the LSP, it connects the period of variation of the star with its SED in the same manner 
as fundamental and overtone stellar pulsation. 



3.3. RGB P-L relation 

Figure [7] shows how the slopes of the RGB P-L relations on sequences 3, 4, E and D relate 
to those determined for the various classes of AGB candidates. The slope of the RGB-star P-L 
relation is relatively wavelength-independent (qualitatively similar to the behavior seen in 0-rich 
candidates) and is consistently less steep than the AGB-star slope at all wavelengths. This effect 



has been noted by other authors (Ita et al. 2004 Glass et al. 2009), and the persistence of this 



trend at longer wavelengths lends support to the conclusion that the populations above and below 
the TRGB on these sequences are indeed dominated by different classes of star. 



Stellar evolutionary models (e.g. Castellani et al. 2003 ) predict (for a limited mass range) an 



offset in variability period at a given luminosity between AGB stars and those on the first ascent 
of the RGB of 6logP ss 0.03, with RGB stars having the longer period. Kiss & Bedding (2003) 



observed this effect in their OGLE based survey, and we confirm it using MACHO. This effect is 
not immediately visually apparent in our sample due to crowding (Figure [2]). However, by solving 
our liner P-L relationships for logP at a fixed [3.6] = 12.0, we find 6logP = 0.05 between the stars 
we define as RGB candidates and those classified as AGBs on sequences 3, 4, and D. 



Nicholls et al. (2010) examined the orbital parameters of a sample of sequence E binary stars, 



and demonstrated that their sample did indeed consist of ellipsoidal binary systems, in contrast to a 
sample of field red giants and a sample of stars from sequence D ( Nicholls et al.||2009 ). In addition, 
they looked for, but did not find, evidence of an IR excess in the colors of sequence E stars, an 
indication of mass loss. Figure [7| shows that the P-L relations of the RGB stars on Sequence E and 
D behave, as a function of wavelength, exactly like the the RGB stars on sequences 3 and 4. The 
lack of change in the slope of the P-L relation indicates that period is not coupled to mass loss in 
these stars. 

The P-L relations of RGB stars on sequences E & D are systematically less steep than those 
of the RGB stars on the pulsational sequences 3 and 4. The slope may be less steep because of 
the photometric definition of an RGB star used in this study, see Appendix [Al Nevertheless, we 
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confirm the null result of Nicholls et al. ( 2010 ) in that we see no evidence of any difference in the 



RGB P-L relations as a function of wavelength for sequence E & D compared to those of sequence 
3 and 4. 



4. DISCUSSION 
4.1. Comparison to Previous Work 

The period-luminosity relationship for AGB stars has been examined by many authors over 



the past several decades (e.g. Feast et al. 1989; Hughes & Wood 1990 Wood et al. 1999 Glass et 
ar]|2009 ) . Previous studies have often relied on a few tens to a few hundreds of evolved LPVs with 
brightnesses often measured in the Kg band. The present work is notable for having the largest 
sample of LPVs with photometric data measured in the most IR bands to date. This has allowed 
us to reduce the uncertainties in these relations by a factor of ~5. 



4-1.1. Comparison to Feast et al. (1989) 



Feast et al. (1989) observed 49 evolved LPVs in the LMC and constructed P-L relations in 



the J, H and K^ bands, as well as bolometric magnitudes, which were used in several subsequent 
studies (e.g. Vassiliadis fc Woo"dl|1993 Whitelock et al.pOOS ). In the K^ band, they find a slope of 
—3.47 lb 0.19 fits the 0-rich stars in their sample, consistent with our slope of —3.31 it 0.04. Our 
fit to the C-rich AGB candidates in the K^ band is also consistent with their derived relation. We 
find a slope of —3.16 it 0.04 compared to their value of —3.30 it 0.40. They do not consider their 
fits to C-rich Miras in the J and H bands to be useful. 



4.1.2. Comparison to Ita et al. (2004) 



Ita et al. (2004) used variability information from the OGLE survey of the LMC (Udalski et 



ar]|1997 ) to assemble a dataset of 35 000 variable stars in the LMC. This sample size allowed them 
to confirm the existence of a distinct population of variable stars below the TRGB, concentrated 



on sequences 3 and 4 (their sequences A and B) first noted by Kiss & Bedding (2003). They reserve 



judgment on the interpretation, put forth by Kiss & Bedding (2003) and used in this paper, that 



this population consists of a large number of first-ascent RGB stars with some AGB contamination 
Most of the Xg P-L relations fits derived in our study (Table pi) are consistent with those in Ita 



et al. (2004). They obtain a slope of of —3.369 it 0.099 for the C-rich stars on sequence 1, consistent 



with our slope (§ 4.1.1). Our most discrepant fits are those to stars brighter than the TRGB in 
sequences 3 and 4. We find a slope of —3.73 it 0.04 for the 1434 0-rich AGB stars we identify 



on sequence 4. Ita et al. (2004) do not differentiate between 0-rich and C-rich AGB stars here. 
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but they find a slope of —3.289 it 0.047 to all 510 stars brighter than the TRGB on this sequence. 
Similarly, on sequence 3 our 0-rich AGB stars are best fit with a slope of —3.85 it 0.04 compared 

This is most likely due to slight differences in sequence 



(2004) 



to -3.356 ± 0.052 in |Ita et al 

definition. Sequence 3 is almost blended with sequence 2, and precisely where the dividing line 

between these two sequences is drawn can greatly impact the derived fit. 



Our study provides a valuable complement to that of Ita et al. (2004), due to the fact that our 



sample is of comparable size but is derived from completely independent observations. 



4-1.3. Comparison to Glass et al. (2009) 



Glass et al. (2009) did their own reduction of a subset of the SAGE and MACHO datasets 



to compute P-L relationships for AGB stars in multiple bands in the IR. Our current work differs 
most noticeably in sample size, their sample consisting of ~1800 stars and ours ~30 000. Precise 
comparison of our derived P-L relations is complicated, however, by several differences in how we 



categorize our datasets. The sample of Glass et al. (2009) is not large enough to distinguish what 



we identify as sequences 2 &: 3, combining them into a single sequence B. See Table [3] for the 



different naming conventions used. Finally, Glass et al. (2009) only fit P-L relations to stars they 
classify as 0-rich, defined as a star with (J — K^) < 1.6, a redder definition than our own (detailed 
2.2). The reddest star that we classify as 0-rich has (J — K^) = 1.4. 



in section 



With these differences in mind, our results show reasonable agreement with those of Glass et 



al. (2009). On sequence 1 for example, our fits for the AGB stars are consistent within the errors 



in nearly every band from Ks to 8 /im. The notable exception to this broad consistency is in the 



IRAC [5.8] band on sequence 1, which we discuss above (§ 3.1) as having an anomalously less steep 



P-L relation. This effect is not observed in Glass et al. ( |2009 ), a result we attribute to sample size. 
Their fit is based on 41 AGB stars on sequence 1, whereas ours is based on 1812. 



4.1.4. Comparison to Vijh et al. (2009) 



By comparing epoch 1 and epoch 2 fiuxes from the SAGE point source catalog, Vijh et al. 



(2009) identified 1967 variable sources in the LMC. Having only two data points, they were not 



able to determine specific variability parameters such as period or amplitude for these sources. We 
identify 731 of these sources in the SAGE single frame mosaic photometry archive with MACHO 



detections. Of the remaining sources, with no MACHO counterpart, Vijh et al. (2009) identified 



41 as 0-rich, 60 as C-rich, and 540 as extreme AGB stars. Only 1.5% of our sample is classified 



as extreme. On the other hand, 84% of the AGB variables identified by Vijh et al. (2009) without 



a MACHO counterpart are classified as extreme AGB stars. By requiring our sources to have a 
MACHO detection, this effect is expected. The MACHO survey focused on the bar of the LMC, and 



extreme AGB candidates are not as concentrated in the bar as are O- and C-rich AGB stars (Blum 



IT 



et al. 2006, fig. 1). Furthermore, because the MACHO survey used two optical bandpasses, we 



expect our sample to miss the reddest, most enshrouded evolved stars in the SAGE catalog, which 



are not visible in the optical. The variables identified by Vijh et al. (2009) are indeed systematically 



redder than the sample identified here. Figure [9] presents a CMD of all the AGB candidate sources 



in this paper, using the SAGE [3.6] and [8.0] bands. The variables detected by Vijh et al. (2009) 



are overlaid as large diamonds. Figure 10 presents a histogram comparing the [3.6]-[8.0] colors of 



the extreme sources in our study with those of Vijh et al. (2009). 



4.1-5. Summary of Com,parison to Other Surveys 

Overall, our quantitative results are consistent with previous studies of the P-L relation in 
AGB stars. Differences that exist are small, and can be attributed to differences in sample size. 
Our current work exhibits two primary strengths, sample size and wavelength coverage. Consisting 
of a sample of 30 000 evolved stars with magnitudes measured in 8 different near-IR bands, our 



sample is far larger than most previous studies. Ita et al. (2004) use a similar sized dataset, but 



lack the wavelength coverage the SAGE survey gives this work. Similarly, Glass et al. ( 2009 ) also 



use the SAGE catalog to investigate the dependence of the P-L relation on wavelength, but use 
only a small fraction of the full SAGE source list. We combine the individual strengths of these 
two works to produce a comprehensive view of evolved star variability in the IR. 



Vijh et al. (2009) also utilized the entire SAGE catalog to probe variability. We expand on this 



work with the inclusion of the MACHO catalog, which allows us to precisely determine the period 
and amplitude of the variation of these stars. However, this additional information introduces 
two selection biases to our sample that do not affect the Vijh sample. The MACHO survey was 
conducted in two non-standard optical filters. By requiring that all of our sources have valid 
MACHO detections, we miss the reddest of the evolved stars. In addition, the MACHO survey 
focused on the bar of the LMC, and this spatial bias also affects our sample. 



4.2. Comparison to Models of AGB Evolution 

Figure ^ reproduces the [3.6] vs. period plot of Figure [ll but this time color coded according 
to RGB, 0-rich, C-rich, or extreme AGB classification. By grouping the sources according to this 



classification, we see trends broadly consistent with synthetic AGB evolution codes (e.g. Vassiliadis 



& Wood 1993), which predict that the structural and chemical changes during the star's evolution 



along the AGB should be accompanied by a general trend toward slightly higher luminosities and 
longer pulsational periods. 

Variable stars on their first ascent of the RGB are dimmer than the AGB stars at later stages of 
evolution, and highly concentrated on sequences 3 and 4, the shortest period sequences. Oxygen-rich 
AGB candidates are distributed across all 4 pulsational sequences. AGB candidates photometrically 
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classified as carbon-rich are also found on all sequences, but are much more concentrated on the 
long period sequences (only 29 are on sequence 4, compared to ~2000 on sequence 1). Extreme 
AGB stars, which are heavily enshrouded by circumstellar dust, are found in significant numbers 
only on sequence 1, the longest period sequence. Table [4] details the number of each sub-class of 
star on each sequence. We see clear trends toward brighter magnitudes and longer periods as one 
moves from RGB candidates to 0-rich, C-rich, and finally extreme AGB candidates. Figures 11 



&; 12 illustrate the distributions of our sample in period and [3.6] magnitude. Previous studies 
(e.g. Feast et ar]|1989 Whitelock et al.|[2003 ), have noticed a dramatic steepening of the AGB P-L 
relation at ~420 days. Figure 12 shows that this is precisely the period at which extreme AGB 



candidates come to dominate over all other classes of stars. Note that in figure 12 the extreme 
AGB candidate graph has been magnified to be visible compared to the much larger 0-rich and 
C-rich AGB populations. Extreme AGB candidates do dominate at periods greater than 420 days. 



3.1 , we find that this class of star 



but the effect is smaller than figure 12 shows. As we discuss in 
does obey a steeper P-L relation. 

The slope of the P-L relation is a measure of the brightness contrast between stars with long 
periods and those of the same type with shorter periods. If the P-L relation becomes increasingly 
steeper with longer wavelength, it is an indication that longer-period stars are emitting more of their 
energy at those wavelengths, an indication of a greater amount of circumstellar matter re- processing 



the star's light to the red. More circumstellar matter implies a higher rate of mass loss ( Srinivasan 



et al.|[20()9 and references therein). Thus, we interpret the lack of a relationship between P-L slope 
and wavelength amongst RGB stars and 0-rich AGB stars to indicate that pulsation period does 
not strongly impact mass-loss rates at these stages of stellar evolution. However, amongst heavily 
enshrouded extreme AGB stars, there is a stronger coupling between the two. 



Other authors (e.g. Schultheis et al. 2004 Glass et al. 2009) have examined plots of color vs. 



period in order to ascertain a period at which significant reddening indicates the onset of mass 
loss. Our determination of the slope of the P-L relation at multiple wavelengths contains and 
quantifies these other studies. Quantitative period-color relations can be easily derived from the 
fits given in this paper. Given P-L fits in two bands, mi and m-2, we obtain the color-period relation 
{mi —1712) = {Al — A2) *logF + {Bi — B2) with A and B defined as in Eqn. 3.1 Groenewegen (2006) 
suggests the [3.6] — [4.5] color as a good indicator of the mass-loss rate from AGB stars. Through 
the period-color relations that can be derived from this work, such proxies can be readily connected 
to variability parameters, and the results compared to the predictions of theoretical models of AGB 
evolution. We intend to investigate the connections between variability and mass-loss rate using 
the quantitative predictions of such models in future papers. 



5. Conclusions 



We present the largest multi- wavelength IR investigation of the period-luminosity relationship 
of evolved stars in the LMC to date. We find that the slope of the P-L relation followed by O- 
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rich AGB stars and stars on the first ascent of the RGB is relatively independent of wavelength. 
In contrast, we find a strong dependence of slope of the P-L relation with wavelength amongst 
the most evolved, most enshrouded, extreme AGB stars. C-rich AGB stars, which represent an 
intermediate stage of evolution along the AGB, show a slight dependence of slope on wavelength, 
with a possible correlation between variability period and the strength of the 5 //m CO absorption 
feature. 

We find that for accurate characterization of the IR P-L relationship for AGB stars, the K^ 
band is the best choice. It offers a superior combination of sample size; low scatter due to the fact 
the both 0-rich and C-rich AGB stars obey identical P-L relations when measured in this band; 
and further applicability, in that once the P-L relation is determined in the K^^ band, the slope of 
the relation for 0-rich AGB stars and RGB stars has been determined at all near-IR wavelengths. 
Slopes of P-L relations for evolved stars in the K^ band can be found in Table [6| 

There is a systematic trend towards steeper slopes in the magnitude vs. logP relation as 
one looks at categories of stars at later and later stages of stellar evolution. This trend would 
be consistent with a non-linear P-L relation if all variable evolved stars were examined together, 
without photometrically determined classifications. 
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Table 1. Dataset Populations 



Source O-rich C-rich Extreme AGB stars Total AGB stars RGB stars Total Sources 



iFraser et al.|||2008|'' N/A N/A 

Srinivasan ct al. (2009)*" 17 958 5179 

|Vijhct al.| | |2009| >° 353 426 

This work 12 172 4455 



N/A 


N/A 


N/A 


56 453 


1428 


24 565 


N/A 


24 565 


820 


1599 


N/A 


1967 


432 


17 059 


13 688 


30 747 



Note. — Number of sources used in recent studies of AGBs and variable sources in the LMC. Classification of sources as 
RGB, O-rich, etc. is done photometrically, using eqns. [l]-[3] 

' ^Fraser et al.| | |2008[ l did not classify their AGB candidates as O-rich, C-rich, etc. 

''Did not use RGB classification 
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Table 3. Period-Luminosity Sequence Conventions 



iFraser et al.lj2005|l 



Sequence 1 
Sequence 2 
Sequence 3 
Sequence 4 
Sequence D 
Sequence E 



Wood et al. 



1 1999 I 



Ita et al. 



(20041 



Glass et al. 



I 2009 I 



Theoretical Explanation 



Sequence C Sequence C Sequence C Pulsation in Fundamental Mode 

Sequence B Sequence C' Sequence B First/Second Overtone Pulsation 

Sequence B Sequence B Sequence B First/Second Overtone Pulsation 

Sequence A Sequence A Sequence A Third Overtone Pulsation 

Sequence D Sequence D Sequence D Undetermined 

Sequence E Sequence E N/A Ellipsoidal binaries 



Note. — Naming conventions used to refer to the various period-luminosity sequences first identified by |Wood et al.| 
l |1999| |. This paper follows the convention of [Eraser et al.| ( |2005^ , because the numeric naming scheme more closely aligns 
with the theoretical explanation. 



Table 4. Sequence Populations 



Sequence 


O-rich 


C-rich 


Extreme 


RGB 


Sequence 1 


2229 


1817 


371 





Sequence 2 


2292 


891 


1 





Sequence 3 


1707 


56 


1 


1992 


Sequence 4 


1434 


29 





2068 


Sequence D 


3364 


971 


18 


5611 


Sequence E 











1727 



Note. — Populations of sources of various types on 
each sequence. Amongst the stellar pulsational sequences, 
RGB stars are concentrated on the short period sequences 
(1-4), O-rich AGB candidates are equally distributed 
across all sequences, and C-rich and extreme AGB candi- 
dates are concentrated on the long period sequences. 
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Table 5. [3.6] vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


r2 


Number 






Sequence 


4 








RGB 


-3.67 ± 0.05 


17.45 ± 0.08 




0.127 


0.687 


2064 


O-rich 


-3.90 ± 0.04 


17.58 ± 0.06 




0.121 


0.886 


1432 


C-rich 


-4.14 ± 0.36 


18.04 ± 0.58 




0.144 


0.829 


29 


All Stars 


-4.56 ± 0.02 


18.68 ± 0.03 
Sequence 


3 


0.138 


0.935 


3525 


RGB 


-3.32 ± 0.05 


17.57 ± 0.08 




0.122 


0.704 


1983 


O-rich 


-4.09 ± 0.04 


18.62 ± 0.07 




0.182 


0.867 


1704 


C-rich 


-3.81 ± 0.27 


18.04 ± 0.51 




0.248 


0.788 


56 


All Stars 


-4.51 ± 0.02 


19.45 ± 0.03 
Sequence 


2 


0.167 


0.939 


3744 


O-rich 


-4.24 ± 0.02 


19.51 ± 0.04 




0.153 


0.942 


2289 


C-rich 


-4.25 ± 0.05 


19.32 ± 0.11 




0.189 


0.890 


887 


All Stars 


-4.54 ± 0.02 


20.04 ± 0.03 
Sequence 


1 


0.177 


0.960 


3176 


O-rich 


-3.41 ± 0.04 


18.89 ± 0.10 




0.271 


0.729 


2221 


C-rich 


-3.77 ± 0.05 


19.35 ± 0.12 




0.251 


0.777 


1813 


Extreme AGB 


-4.27 ± 0.19 


20.37 ± 0.49 




0.336 


0.586 


371 


All Stars 


-4.22 ± 0.02 


20.55 ± 0.06 
Sequence 


D 


0.304 


0.875 


4405 


RGB 


-2.68 ± 0.03 


19.33 ± 0.07 




0.211 


0.608 


5588 


O-rich 


-3.58 ± 0.04 


21.39 ± 0.10 




0.217 


0.752 


3354 


C-rich 


-3.98 ± 0.08 


22.17 ± 0.25 




0.299 


0.701 


969 


Extreme AGB 


-4.71 ± 1.40 


23.72 ± 4.28 




0.545 


0.414 


18 


All Stars 


-4.16 ± 0.02 


23.12 ± 0.04 
Sequence 


E 


0.288 


0.863 


9929 


RGB 


-2.62 ± 0.02 


18.52 ± 0.05 




0.274 


0.894 


1718 


All Stars 


-2.62 ± 0.02 


18.52 ± 0.05 




0.274 


0.894 


1718 



Note. — Parameters for the linear fit to the IRAC [3.6] log P-magnitudc relation. 
Identical tables for the J, H, [4.5], [5.8], [8.0] and [24] bands are available online. 
"Residual" is the standard error of the residuals to the fit. The r^ column lists 
the correlation coefBcicnt, and "number" lists the number of sources used to derive 
each fit. 
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Table 6. K^ vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


,2 


Number 






Sequence 


4 








RGB 


-3.43 ± 0.05 


17.22 ± 0.07 




0.111 


0.677 


2054 


O-rich 


-3.73 ± 0.04 


17.43 ± 0.06 




0.126 


0.867 


1434 


C-rich 


-3.73 ± 0.32 


17.58 ± 0.52 




0.134 


0.833 


29 


All Stars 


-4.46 ± 0.02 


18.65 ± 0.03 
Sequence 


3 


0.135 


0.929 


3517 


RGB 


-3.12 ± 0.05 


17.36 ± 0.08 




0.117 


0.684 


1981 


O-rich 


-3.85 ± 0.04 


18.33 ± 0.07 




0.182 


0.846 


1707 


C-rich 


-3.61 ± 0.19 


17.99 ± 0.35 




0.182 


0.875 


56 


All Stars 


-4.33 ± 0.02 


19.27 ± 0.03 
Sequence 


2 


0.166 


0.932 


3744 


O-rich 


-4.06 ± 0.02 


19.33 ± 0.04 




0.143 


0.948 


2292 


C-rich 


-3.69 ± 0.04 


18.61 ± 0.09 




0.155 


0.899 


891 


All Stars 


-3.93 ± 0.01 


19.10 ± 0.03 
Sequence 


1 


0.149 


0.963 


3183 


O-rich 


-3.31 ± 0.04 


18.87 ± 0.09 




0.272 


0.760 


2218 


C-rich 


-3.16 ± 0.04 


18.40 ± 0.11 




0.230 


0.744 


1817 


Extreme AGB 


-2.56 ± 0.35 


17.43 ± 0.91 




0.444 


0.147 


312 


All Stars 


-3.34 ± 0.02 


18.90 ± 0.05 
Sequence 


D 


0.293 


0.848 


4347 


RGB 


-2.60 ± 0.03 


19.26 ± 0.07 




0.197 


0.622 


5611 


O-rich 


-3.40 ± 0.03 


21.04 ± 0.09 




0.197 


0.767 


3364 


C-rich 


-3.58 ± 0.06 


21.50 ± 0.18 




0.220 


0.774 


971 


Extreme AGB 


-4.41 ± 0.64 


24.04 ± 1.96 




0.246 


0.748 


18 


All Stars 


-3.81 ± 0.01 


22.34 ± 0.04 
Sequence 


E 


0.236 


0.873 


9964 


RGB 


-2.54 ± 0.02 


18.44 ± 0.05 




0.246 


0.893 


1706 


All Stars 


-2.54 ± 0.02 


18.44 ± 0.05 




0.246 


0.893 


1706 



Note. 



Identical to table [s] with magnitude measured in the Ks band 
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• Sequence 4 

• Sequence 3 

• Sequence 2 

• Sequence 1 

• Sequence D 

• Sequence E 







Period (Days) 



Fig. 1. — Period-luminosity sequences for evolved stars in the LMC. Sequence 1 consists of stars 
pulsating in the fundamental mode, while sequences 2-4 are higher order pulsational modes. Se- 
quence E consists of ellipsoidal binary systems, and the mechanism responsible for the variation on 



sequence D is not known. The naming convention follows that of Fraser et al. (2008). This figure 
is best appreciated in color, available in the online version. 
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• Oxygen-Rich AGB 

• Carbon-Rich AGB 

• Extreme AGB 







10 



100 
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1000 



Fig. 2. — A reproduction of Figure [T| but color coded according to the photometrically determined 
classification of the source. See § |2.2| for details. This figure is best appreciated in color, available 
in the online version. 
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600 



400 



200 




11 10 

[3.6] magnitude 



Fig. 3. — Histogram of [3.6] magnitude, color-coded by period-luminosity sequence. The RGB 
population is clearly visible below 12*'^ magnitude on sequences 3, 4, E and D, becoming less 
distinct on sequences 1 and 2. Similarly, C-rich AGB stars, brighter than magnitude 10.5, are 
much more prominent on sequences 1 and 2. 
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Fig. 4. — Four versions of the period luminosity diagram of figure [2j using different bands as the 
luminosity proxy. Top row, L to R: the 2MASS K^ band and the Spitzer IRAC 4.5 /im band from 
the SAGE survey. Bottom row, L to R: the IRAC 8.0 /im band and the MIPS 24 /xm band. The 
O- and C-rich AGB candidates maintain the same relative position in all the bands. The extreme 
sources do not look particularly bright in the K^ band, but further into the IR, are revealed as the 
brightest sources in our sample. The RGB stars are almost entirely absent from the MIPS 24 ^um 
diagram. 
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8- 






10 = 



12- 



14- 



I III 



Oxygen-Rich AGB 
Carbon-Rich AGB 
Extreme AGB 




100 



1000 



Period (Days) 



Fig. 5. — Linear least-squares fits to the period-[3.6] magnitude relation for the sources assigned 
to sequence 1, the fundamental mode pulsators. Fits were done to each sub-group individually. 
Sequence 1 is the only sequence with a significant number of extreme AGB sources on it. Notice 
the smooth progression from shallow to steep slopes for the 3 different classes of AGB stars. 



27 



Q. 



CO 




pel [4.5] 

Wavelength 



Fig. 6. — The slopes of the period-magnitude relations for each type of star in sequences 1 (top) 
and 2 (bottom), plotted for all the SAGE IR wavelengths. Steeper slopes are more negative, and 
are towards the bottom of this plot. A steeper slope is an indication that longer period variables 
have brightened relative to the shorter period stars on the same sequence. 
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Fig. 7. — Similar to figure [6j but showing the 4 P-L sequences not shown there. Top row: Left: 
Sequence 3 Right: Sequence 4. Bottom row: Left: Sequence D Right: Sequence E 
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Fig. 8. — Median SEDs for each of the four population classes usetd in this paper, first ascent RGB 
(top left), 0-rich AGB candidates (top right), C-rich AGB candidates (bottom left) and extreme 
AGB candidates (bottom right). Error bars represent the spread in the population, not measure- 
ment uncertainty. The utilization of longer wavelength IRAC bands allows us to more accurately 
measure the intrinsic brightness of these sources, particularly the extreme AGB candidates, whose 
SEDs peak in the 3-8 /xm range. 
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Fig. 9.— [8.0] vs. [3.6] - [8.0] CMD of the AGB candidates detected by [Vijh et al] ( |2009D overlaid 
atop the CMD for the sources used in this paper, detected by both SAGE and MACHO. The 



sources from Vijh et al. (2009) appear as large, open diamonds. Note the long tail of very red 
extreme AGB candidates ([3.6] - [8.0]) > 4 detected by SAGE but not MACHO. 
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Fig. 10. — Histogram of [3.6] — [8.0] color of the sources classified as extreme AGB candidates in 
this work and in Vijh et al. ( |2009 ). Because we require a valid MACHO (optical) detection, we 



miss the reddest and most enshrouded of the Vijh et al. (2009) extreme sources 
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Fig. 11. — Histogram of IRAC 3.6 fim. magnitude for the four types of sources in our sample. The 
extreme AGB sample has been artificially repUcated 5 times to make it visible on the same scale 
as the others. 
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Fig. 12. — Histogram of MACHO primary period for the four types of sources in our sample. The 
extreme AGB sample has been artificially replicated 5 times to make it visible on the same scale 
as the others. 
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A. Period-Luminosity Sequence Definitions 

The stars in our sample are assigned to one of 6 P-L sequences using cuts in i^s-Period space, 
witli a small subset also utilizing J—K^ color infonnation. Only the primary period of these multi- 
periodic objects is used to establish sequence membership. The definitions for the sequences are 
shown graphically in figure 13 and described mathematically in Table [7j Our definitions are based 



on those used by Eraser et al. (2008), but are more restrictive (especially in sequence D, which 
exhibits a very large scatter about the P-L relation) to better trace the population density "cores" 
visible to the eye. This emphasis sacrifices some sample size, but we note that the scatter in our 
P-L relations (of order ~0.15 mag) is dominated by the intrinsic variability of the sources and 
would not be significantly reduced by introducing more sources by using the secondary periods to 
assign the un-sequenced sources to one of the existing P-L sequences. 

We note that ultimately, our P-L sequence definitions are a judgment call, and slightly different 
definitons used by other authors are equally valid and can affect the measured P-L fits. In particular, 
the TRGB cut at a constant Ks = 12.0 leads to a sharp parallelogram shape at the brightest end 
of the RGB population. These corners cause the RGB and 0-rich AGB P-L fits to be less steep 
than the fit to all the stars in a sequence, regardless of class, which is not affected by this selection 
bias (Table [5]). 

Beginning with these definitions of the P-L sequences, careful study of the panels of Figure [4] 
led us to introduce some additional modifications to properly classify a few anomalous populations. 

A group of 100 C-rich AGB candidates, classified as being part of sequence 1 under criterion 
1 in table [TJ was found to remain sub-luminous to the rest of sequence 1 at all wavelengths longer 
than ifg as well. They are identified and eliminated by removing C-rich AGB candidates from 
sequence 1 with periods > 300 days and Kg > —3.9 * logP + 20.17. 

Finally, we identified a population of 60 extreme AGB candidates so dust-enshrouded that at 
Ks they lay along the main body of sequence D, but at 3.6 /im and longer they were clearly part 
of sequence 1. They were properly classified by adding all extreme AGB candidates assigned to 
sequence D or to no sequence with periods between 250 and 700 days and IRAC [3.6] < 10.5 to 
sequence 1. 

Note that both of these populations also require their members to be photometrically classified 



as either C-rich or extreme AGB candidates, following the definitions detailed in § 2.2 
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Period (Days) 



Fig. 13. — Kg band P-L diagram of our sample used in this study, showing the boundaries which 
define the 6 P-L sequences identified in this work. The populations of these sources are detailed 
in Table |4] The boundaries are described quantitatively in table [7j The blue region in between 
sequences 1 and D categorizes stars with J - Kg > 1.4 as being on sequence 1. 
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B. Mid-IR Period-Magnitude Relations 
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Table 8. J vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


r2 


Number 






Sequence 


4 








RGB 


-3.00 ± 0.06 


17.65 ± 0.08 




0.126 


0.558 


2050 


O-rich 


-3.46 ± 0.04 


18.08 ± 0.07 




0.142 


0.814 


1434 


C-rich 


-3.74 ± 0.33 


18.95 ± 0.53 




0.128 


0.844 


26 


All Stars 


-4.16 ± 0.02 


19.27 ± 0.04 
Sequence 


3 


0.153 


0.899 


3510 


RGB 


-2.60 ± 0.05 


17.56 ± 0.08 




0.125 


0.579 


1979 


O-rich 


-3.51 ± 0.04 


18.82 ± 0.08 




0.195 


0.797 


1707 


C-rich 


-3.31 ± 0.22 


18.79 ± 0.42 




0.214 


0.823 


51 


All Stars 


-3.91 ± 0.02 


19.62 ± 0.03 
Sequence 


2 


0.181 


0.906 


3737 


O-rich 


-3.57 ± 0.02 


19.51 ± 0.04 




0.154 


0.927 


2290 


C-rich 


-2.76 ± 0.06 


18.15 ± 0.12 




0.220 


0.719 


887 


All Stars 


-2.96 ± 0.02 


18.43 ± 0.04 
Sequence 


1 


0.223 


0.887 


3177 


O-rich 


-2.87 ± 0.04 


18.99 ± 0.08 




0.252 


0.747 


2219 


C-rich 


-2.10 ± 0.05 


17.46 ± 0.13 




0.307 


0.460 


1812 


Extreme AGB 


-1.47 ± 0.50 


17.19 ± 1.32 




0.767 


0.027 


311 


All Stars 


-1.95 ± 0.03 


17.06 ± 0.06 
Sequence 


D 


0.500 


0.539 


4342 


RGB 


-2.04 ± 0.03 


18.81 ± 0.07 




0.196 


0.515 


5590 


O-rich 


-3.13 ± 0.03 


21.40 ± 0.09 




0.190 


0.749 


3364 


C-rich 


-2.99 ± 0.06 


21.25 ± 0.19 




0.240 


0.700 


965 


Extreme AGB 


-2.11 ± 1.81 


19.40 ± 5.55 




0.600 


0.088 


16 


All Stars 


-3.15 ± 0.01 


21.64 ± 0.04 
Sequence 


E 


0.248 


0.828 


9935 


RGB 


-2.24 ± 0.02 


18.73 ± 0.04 




0.236 


0.878 


1703 


All Stars 


-2.24 ± 0.02 


18.73 ± 0.04 




0.236 


0.878 


1703 



Note. 



Identical to table [s] with magnitude measured in the J band 
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Table 9. H vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


r2 


Number 






Sequence 


4 








RGB 


-3.23 ± 0.05 


17.14 ± 0.07 




0.112 


0.648 


2048 


O-rich 


-3.55 ± 0.04 


17.37 ± 0.07 




0.132 


0.842 


1432 


C-rich 


-3.37 ± 0.33 


17.32 ± 0.53 




0.145 


0.798 


28 


All Stars 


-4.26 ± 0.02 


18.56 ± 0.03 
Sequence 


3 


0.139 


0.919 


3508 


RGB 


-2.92 ± 0.05 


17.24 ± 0.07 




0.115 


0.668 


1975 


O-rich 


-3.61 ± 0.04 


18.14 ± 0.08 




0.188 


0.818 


1704 


C-rich 


-3.36 ± 0.19 


17.90 ± 0.36 




0.190 


0.863 


52 


All Stars 


-4.06 ± 0.02 


19.02 ± 0.03 
Sequence 


2 


0.169 


0.923 


3731 


O-rich 


-3.75 ± 0.02 


19.00 ± 0.04 




0.145 


0.941 


2291 


C-rich 


-3.10 ± 0.05 


17.83 ± 0.11 




0.187 


0.814 


888 


All Stars 


-3.35 ± 0.02 


18.28 ± 0.03 
Sequence 


1 


0.181 


0.935 


3179 


O-rich 


-3.06 ± 0.04 


18.56 ± 0.08 




0.267 


0.753 


2218 


C-rich 


-2.52 ± 0.05 


17.40 ± 0.12 




0.268 


0.605 


1814 


Extreme AGB 


-1.70 ± 0.43 


16.30 ± 1.12 




0.620 


0.048 


314 


All Stars 


-2.53 ± 0.02 


17.45 ± 0.06 
Sequence 


D 


0.398 


0.720 


4346 


RGB 


-2.41 ± 0.03 


18.95 ± 0.07 




0.194 


0.596 


5590 


O-rich 


-3.19 ± 0.03 


20.72 ± 0.09 




0.189 


0.757 


3361 


C-rich 


-3.19 ± 0.06 


20.83 ± 0.18 




0.220 


0.745 


967 


Extreme AGB 


-2.30 ± 1.07 


18.57 ± 3.27 




0.375 


0.249 


16 


All Stars 


-3.46 ± 0.01 


21.60 ± 0.04 
Sequence 


E 


0.231 


0.860 


9934 


RGB 


-2.46 ± 0.02 


18.43 ± 0.04 




0.237 


0.896 


1700 


All Stars 


-2.46 ± 0.02 


18.43 ± 0.04 




0.237 


0.896 


1700 



Note. 



Identical to table [s] with magnitude measured in the H band 
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Table 10. [4.5] vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


r2 


Number 






Sequence 


4 








RGB 


-3.39 ± 0.05 


17.18 ± 0.07 




0.110 


0.711 


2065 


O-ricii 


-3.86 ± 0.04 


17.66 ± 0.06 




0.121 


0.885 


1434 


C-rich 


-4.14 ± 0.35 


18.17 ± 0.57 




0.142 


0.836 


29 


All Stars 


-4.42 ± 0.02 


18.61 ± 0.03 
Sequence 


3 


0.128 


0.938 


3528 


RGB 


-3.04 ± 0.04 


17.23 ± 0.07 




0.104 


0.728 


1992 


O-rich 


-4.09 ± 0.04 


18.77 ± 0.07 




0.172 


0.877 


1706 


C-rich 


-3.87 ± 0.27 


18.26 ± 0.51 




0.246 


0.796 


56 


All Stars 


-4.36 ± 0.02 


19.31 ± 0.03 
Sequence 


2 


0.155 


0.943 


3755 


O-rich 


-4.16 ± 0.02 


19.43 ± 0.04 




0.154 


0.936 


2290 


C-rich 


-4.05 ± 0.05 


18.99 ± 0.12 




0.204 


0.862 


891 


All Stars 


-4.44 ± 0.02 


19.93 ± 0.03 
Sequence 


1 


0.180 


0.955 


3181 


O-rich 


-3.35 ± 0.05 


18.80 ± 0.10 




0.270 


0.708 


2227 


C-rich 


-3.56 ± 0.05 


18.96 ± 0.12 




0.265 


0.750 


1816 


Extreme AGB 


-4.87 ± 0.22 


21.55 ± 0.58 




0.402 


0.560 


371 


All Stars 


-4.16 ± 0.03 


20.47 ± 0.06 
Sequence 


D 


0.367 


0.854 


4414 


RGB 


-2.41 ± 0.03 


18.72 ± 0.07 




0.214 


0.562 


5610 


O-rich 


-3.45 ± 0.04 


21.14 ± 0.11 




0.242 


0.703 


3361 


C-rich 


-3.87 ± 0.09 


21.93 ± 0.27 




0.338 


0.646 


969 


Extreme AGB 


-5.24 ± 1.95 


25.04 ± 5.94 




0.690 


0.312 


18 


All Stars 


-4.01 ± 0.02 


22.82 ± 0.05 
Sequence 


E 


0.316 


0.838 


9958 


RGB 


-2.56 ± 0.02 


18.46 ± 0.05 




0.275 


0.891 


1727 


All Stars 


-2.56 ± 0.02 


18.46 ± 0.05 




0.275 


0.891 


1727 



Note. 



Identical to tablep] with magnitude measured in the [4.5] band 
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Table 11. [5.8] vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


r2 


Number 






Sequence 


4 








RGB 


-3.48 ± 0.05 


17.19 ± 0.07 




0.105 


0.712 


2063 


O-ricii 


-3.90 ± 0.04 


17.61 ± 0.06 




0.118 


0.891 


1432 


C-rich 


-4.27 ± 0.34 


18.25 ± 0.54 




0.134 


0.858 


29 


All Stars 


-4.52 ± 0.02 


18.63 ± 0.03 
Sequence 


3 


0.126 


0.942 


3524 


RGB 


-3.22 ± 0.04 


17.40 ± 0.07 




0.112 


0.723 


1985 


O-rich 


-4.20 ± 0.04 


18.83 ± 0.07 




0.170 


0.884 


1705 


C-rich 


-3.88 ± 0.27 


18.11 ± 0.51 




0.248 


0.794 


56 


All Stars 


-4.53 ± 0.02 


19.47 ± 0.03 
Sequence 


2 


0.158 


0.946 


3747 


O-rich 


-4.30 ± 0.02 


19.56 ± 0.04 




0.160 


0.935 


2288 


C-rich 


-3.72 ± 0.06 


18.16 ± 0.14 




0.234 


0.794 


887 


All Stars 


-4.45 ± 0.02 


19.80 ± 0.04 
Sequence 


1 


0.190 


0.948 


3175 


O-rich 


-3.41 ± 0.05 


18.77 ± 0.10 




0.287 


0.697 


2223 


C-rich 


-3.22 ± 0.05 


18.04 ± 0.13 




0.308 


0.664 


1812 


Extreme AGB 


-5.35 ± 0.30 


22.44 ± 0.77 




0.509 


0.470 


369 


All Stars 


-4.02 ± 0.03 


20.04 ± 0.07 
Sequence 


D 


0.439 


0.820 


4404 


RGB 


-2.54 ± 0.03 


18.94 ± 0.08 




0.230 


0.557 


5594 


O-rich 


-3.49 ± 0.04 


21.13 ± 0.12 




0.259 


0.690 


3355 


C-rich 


-3.70 ± 0.10 


21.27 ± 0.29 




0.374 


0.589 


970 


Extreme AGB 


-5.43 ± 2.30 


25.23 ± 7.02 




0.781 


0.258 


18 


All Stars 


-4.20 ± 0.02 


23.18 ± 0.05 
Sequence 


E 


0.333 


0.841 


9937 


RGB 


-2.63 ± 0.02 


18.52 ± 0.05 




0.292 


0.877 


1711 


All Stars 


-2.63 ± 0.02 


18.52 ± 0.05 




0.292 


0.877 


1711 



Note. 



Identical to tablep] with magnitude measured in the [5.8] band 
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Table 12. [8.0] vs. logP linear fit parameters 



Classification 


Slope 


Intercept 


Residual 


r2 


Number 






Sequence 


4 








RGB 


-3.56 ± 0.05 


17.29 ± 0.08 




0.122 


0.682 


2044 


O-ricii 


-4.04 ± 0.04 


17.79 ± 0.06 




0.123 


0.892 


1429 


C-rich 


-4.50 ± 0.33 


18.58 ± 0.54 




0.140 


0.871 


29 


All Stars 


-4.62 ± 0.02 


18.76 ± 0.03 
Sequence 


3 


0.136 


0.940 


3502 


RGB 


-3.23 ± 0.05 


17.40 ± 0.08 




0.131 


0.684 


1978 


O-rich 


-4.43 ± 0.04 


19.19 ± 0.07 




0.179 


0.885 


1704 


C-rich 


-4.12 ± 0.34 


18.39 ± 0.65 




0.324 


0.734 


55 


All Stars 


-4.68 ± 0.02 


19.69 ± 0.03 
Sequence 


2 


0.174 


0.942 


3738 


O-rich 


-4.49 ± 0.03 


19.83 ± 0.05 




0.198 


0.911 


2284 


C-rich 


-4.58 ± 0.08 


19.64 ± 0.17 




0.275 


0.800 


889 


All Stars 


-5.05 ± 0.02 


20.81 ± 0.05 
Sequence 


1 


0.250 


0.933 


3173 


O-rich 


-3.45 ± 0.06 


18.74 ± 0.13 




0.344 


0.612 


2221 


C-rich 


-3.95 ± 0.06 


19.40 ± 0.14 




0.317 


0.716 


1816 


Extreme AGB 


-5.77 ± 0.33 


23.09 ± 0.86 




0.546 


0.451 


371 


All Stars 


-4.68 ± 0.03 


21.28 ± 0.08 
Sequence 


D 


0.473 


0.812 


4408 


RGB 


-2.47 ± 0.03 


18.74 ± 0.09 




0.265 


0.479 


5561 


O-rich 


-3.55 ± 0.05 


21.19 ± 0.14 




0.331 


0.605 


3361 


C-rich 


-4.27 ± 0.12 


22.61 ± 0.35 




0.433 


0.575 


969 


Extreme AGB 


-5.74 ± 2.44 


25.72 ± 7.43 




0.836 


0.257 


18 


All Stars 


-4.50 ± 0.02 


23.93 ± 0.06 
Sequence 


E 


0.411 


0.808 


9909 


RGB 


-2.55 ± 0.03 


18.29 ± 0.07 




0.349 


0.787 


1631 


All Stars 


-2.55 ± 0.03 


18.29 ± 0.07 




0.349 


0.787 


1631 



Note. 



Identical to tablep] with magnitude measured in the [8.0] band 
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Table 13. [24] vs. logP linear fit parameters 



Classification 


Slope 


Intercept 




Residual 


2 

r 


Number 






Sequence 


4 








RGB 


-3.21 ± 2.25 


14.72 ± 3.24 




0.591 


0.061 


33 


O-rich 


-0.79 ± 0.27 


11.57 ± 0.45 




0.344 


0.048 


177 


All Stars 


-0.25 ± 0.23 


10.58 ± 0.39 
Sequence 


3 


0.416 


0.005 


213 


RGB 


3.71 ± 1.93 


4.15 ± 3.11 




0.666 


0.091 


39 


O-rich 


-3.02 ± 0.17 


15.82 ± 0.33 




0.429 


0.401 


471 


C-rich 


-1.38 ± 0.77 


12.20 ± 1.52 




0.449 


0.113 


27 


All Stars 


-1.97 ± 0.14 


13.69 ± 0.28 
Sequence 


2 


0.497 


0.257 


538 


O-rich 


-1.45 ± 0.18 


12.75 ± 0.37 




0.537 


0.079 


731 


C-rich 


-3.45 ± 0.18 


16.78 ± 0.39 




0.448 


0.328 


758 


AU Stars 


-2.84 ± 0.11 


15.49 ± 0.22 
Sequence 


1 


0.505 


0.328 


1489 


O-rich 


-0.61 ± 0.23 


11.16 ± 0.54 




0.701 


0.008 


800 


C-rich 


-3.38 ± 0.10 


17.67 ± 0.25 




0.416 


0.410 


1587 


Extreme AGB 


-6.34 ± 0.58 


23.98 ± 1.50 




0.752 


0.249 


366 


AU Stars 


-3.01 ± 0.10 


16.65 ± 0.24 
Sequence 


D 


0.728 


0.254 


2753 


RGB 


0.07 ± 0.71 


9.70 ± 1.85 




0.679 


0.000 


129 


O-rich 


-0.65 ± 0.15 


11.63 ± 0.43 




0.651 


0.016 


1240 


C-rich 


-2.53 ± 0.17 


16.84 ± 0.50 




0.480 


0.227 


789 


Extreme AGB 


-3.04 ± 4.76 


16.39 ± 14.65 




0.990 


0.028 


16 


AU Stars 


-1.43 ± 0.09 


13.75 ± 0.26 
Sequence 


E 


0.649 


0.107 


2174 


RGB 


-1.26 ± 0.71 


11.89 ± 1.58 




1.062 


0.084 


36 


AU Stars 


-1.26 ± 0.71 


11.89 ± 1.58 




1.062 


0.084 


36 



Note. — Identical to table Is] with magnitude measured in the [24] band 



